Thermal fluctuations and longitudinal relaxation of single-domain magnetic particles 

at elevated temperatures 
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We present numerical and analytical results for the swiching times of magnetic nanoparticles with 
uniaxial anisotropy at elevated temperatures, including the vicinity of T c . The consideration is based 
in the Landau-Lifshitz-Bloch equation that includes the relaxation of the magnetization magnitude 
M. The resulting switching times are shorter than those following from the naive Landau-Lifshitz 
equation due to (i) additional barrier lowering because of the reduction of M at the barrier and (ii) 
critical divergence of the damping parameters. 
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The theory of thermal fluctuations of small magnetic 
particles is one of the fundamental issues of modern mi- 
cromagnetics. The conditions at which the particle be- 
comes superparamagnetic define the thermal stability of 
the magnetized state and, therefore, is also valuable for 
technological application such as magnetic recording 
The basis of the theory has been introduced by Brown 
:2| who suggested to include thermal fluctuations into 
the Landau-Lifshitz (LL) dynamical equation as formal 
random fields whose properties are defined by the equilib- 
rium solution of the correspondent Fokker-Planck (FP) 
equation. He also derived the Arrhenius-Neel formula to 
describe the relaxation rate in the axially symmetric case 
of single-domain particles which was lately generalized 
to the presence of external field 0, 0] . Since the paper 
of Chantrell and Lyberatos this Langevin-dynamics 
approach of Brown has been brought to numerical micro- 
magnetics to model the thermal properties of an ensem- 
ble of interacting particles and, more generally, of ferro- 
magnetic materials, interpreting the micromagnetic dis- 
cretization elements as small particles. Generally speak- 
ing, the LL equation used in these simulations is a low- 
temperature approximation only. Recently a generalized 
Landau-Lifshitz-Bloch (LLB) equation for a ferromag- 
net 0, 0] was derived which is valid for all temperatures 
and includes the longitudinal relaxation. The deviations 
of the LLB dynamics from the LL dynamics should be 
pronounced at high temperatures, especially close to the 
Curie temperature T c . The validity of this approach has 
been confirmed by the measurements of the domain-wall 
mobility in crystals of Ba-and Sr-hexaferrites close to T c 



has been suggested pLl| . Therefore, from both funda- 
mental and applied points of view it is necessary to 
consider the micromagnetics of small particles (or mag- 
netic grains) at elevated temperatures. The straightfor- 
ward approach [l2j uses the formalism of the standard 
stochastic LL equation, however with the temperature- 
dependent parameters, i.e., the equilibrium magnetiza- 
tion M e (T) introduced through the mean-field approxi- 
mation (MFA) involving the Brillouin function, and the 
uniaxial anisotropy K{T) through the scaling relation 
K(T) ~ Mg(T). However, this approach becomes in- 
valid at elevated temperatures as it does not incorpo- 
rate the longitudinal relaxational effects. The purpose 
of this Letter is to introduce the theoretical formalism 
of the thermal fluctuations of single-domain particles on 
the basis of the LLB equation which should be valid at 
all temperatures. As a practical example, we consider 
analytically and numerically the thermal switching of a 
FePt particle and discuss the conditions at which the dif- 
ferences between the two formalisms emerge. 

We start with the LLB equation 0, Q augmented by 
the white-noise fields £, £ 1; and C 2 m the form 



n = 7 [n x (H cff + 0] + ^(n- (H cff + d))n 



7a 2 



[nx [nx (H off + C 2 )]], 



(1) 



where n = M/M e (T), is the reduced magnetization, 7 
is the gyromagnetic ratio, a\ and a 2 are dimensionlcss 
longitudinal and transverse damping parameters. The 
effective field H e fl is given by 



Since the proposal of the heat-assisted magnetic 
recording (HAMR) 10] the problem of high-temperature 
thermal magnetization dynamics has become of large 
practical importance. The basic idea of HAMR is to 
write bits of information at elevated temperature (close 
to the Curie temperature, where the switching field is 
small) and store the information at room temperature. 
To achieve a significant areal density advantage, the 
use of high-anisotropy intermetallics such as Llg FePt 



SF 

Hcff = - m = H + H4 + (M e /xn) (1 - n 2 



n, 



(2) 



where F is the free energy density of the single-domain 
particle (cf. Ref. Q), H and are applied and 
anisotropy fields, and \\\ — dM e /dH is the longitudi- 
nal susceptibility. Parameters M e , xp an d &i,2 depend 
on temperature and they are singular at T c . Within the 
MFA-based model, one can use Eq. (4.9) of Ref. with 
K\ = K2 and "f 1 2 => a 1,2, rearranged to the form similar 



2 



to that of of Rcf. 



A 2T 2q 
m e 3T C sinh(2q) 



a 2 = — 



tanh q T 
q 3T C 



(3) 



Here A is a microscopic damping parameter that is tern 
perature dependent but noncritical at T c , 



m e = M e (T)/M e (T = 0) 



(4) 



is the reduced magnetization, and q = 3T c m e /[2(S+l)T). 
One can see that in the vicinity of T c the relaxational 
parameters diverge: a\ = a 2 cc l/M e (T). In accordance 
with this theoretical prediction, ferromagnetic resonance 
measurements on permalloy have shown sharp increase 
of the damping close to the Curie temperature . 

The stochastic fields in the LLB-Langevin equation 
above can be, in fact, introduced in many different ways. 
For instance, one can consider all three fields as uncor- 
related, set some of them to zero, or identify some of 
them with each other. In all these cases one obtains dif- 
ferent LLB-Langevin equations and different stochastic 
trajectories. The physical quantities obtained by averag- 
ing over realizations of £, £ 1; £ 2 are, however, the same 
for all models. The reason is that in all cases one obtains 
the same Fokker-Planck equation (FPE), as was shown 
in Ref. Q for the LL-Langevin equation. The FPE cor- 
responding to Eq. JQ) can be obtained in the same way 
as that for the LL equation (see Appendix of Ref. 0). 
The result has the form of the conservation law 



dt dn 



0, 



(5) 



where / = /(n,t) is the probability density and the prob- 
ability current J reads 



[n x H cff ]/- 



n n- H, 



«-efit 



Ot2 
2 



H, 



T d 
VM e dn 

T t. 



VM e dn 



f 



(6) 



We will illustrate the stochastic dynamics of single- 
domain magnetic particles for the model with the z- 
uniaxial anisotropy, 



F A = (M 2 +M 2 )/(2 X± ), 



(7) 



where x± is the transverse susceptibility that is a con- 
stant within the MFA. We use Eq. JJJ) rather than F A — 
—M 2 /(2x±) to make T c independent of the anisotropy 
and thus to simplify our formalism. Eq. Q could be 
rewritten using a generalization of the micromagnetic 
anisotropy K as F A = {y 2 x + iP\ K (or as F A = —v 2 K), 
where v is the magnetization direction vector, v = M/M 
(see Ref. ||). This is not helpful, however, within the 
approach based on the LLB equation. The problem is 
that K = M 2 /(2x±) is not a constant and even not a 
function of temperature (cf. K = M 2 (T)/(2x±) used 
in Ref. E2)j since the magnetization magnitude M 



can change dynamically during the thermal escape pro- 
cess. It is convenient to scale the free energy density as 
F = {M 2 Jxi) F with F given by^ 



F 



-n • h- 



(8) 



where h = (xj_/M e ) H and a = 2X\\/x±- We also define 
the temperature variable a similarly to Ref. Q : 



VF/T = 2aF, a = VM 2 / {2 Xl T) 



(9) 



We will restrict our consideration to the case H = 0. In 
this case the minima of F are located at n x = n y = 0, 



±1, and Fr, 



and n± = 



= 0. The saddle point of F is n z 
2 = n s , where 



u 



\/l — 7i, a < 1 - __ J (2 - a) /4, a < 1 
0, a>l, ^ sad -\ l/(4a), a>l. 

(10) 

In the limit T -> (i.e., X|| and thus a -> 0) Eq. © 
confines the vector n to the sphere n = |n| = 1, and the 
standard formalism based on the LL equation is recov- 
ered. At nonzero temperatures, a > 0, the magnetization 
changes its magnitude. In our model this effect is maxi- 
mal at the saddle point where the magnetization vector 
is perpendicular to the easy axis. One can visualize the 
trajectory of this vector in the process of thermal activa- 
tion, after averaging over fluctuations, as an ellipsis going 
via the saddle point. At T — T* defined by a — 1 the 
ellipsis degenerates into a line. In the range T* < T < T c 
one has a > 1, and the situation is qualitatively different. 
Here n contracts preserving its direction along the z axis 
and goes through zero at the saddle point, then it grows 
in the opposite direction. These scenarios are very sim- 
ilar to the transformation of the domain wall structure 
with temperature . Obviously the process of thermal 
activation of single-domain magnetic particles cannot be 
described on the basis of the LL equation at elevated 
temperatures. The crucial process of the longitudinal re- 
laxation is captured by Eqs. based on the LLB 
equation. 

The escape rate in the case T <£l AU has the form 



r = Tq exp 



AU 



AU VF 



T 



T 



2aF, 



sad ■ 



ID 



In addition to the dependence a oc M 2 (T) in Eq. © 
that is responsible for the barrier lowering at elevated 
temperatures, there is another source of the barrier low- 
ering described by F sa d in Eq. I |10| l. In particular, the 
value of -F sa d at a = 1 is two times smaller than its low- 
temperature value, a — > 0. The prefactor r in Eq. 111|) 
can be obtained by solving the FPE, Eq. (JSJ, similarly to 
the derivation in Ref. 4] . For a < 1 the result depends on 
a2 only, since in this case the barrier is being crossed by 
the pure rotation of the magnetization vector. For a > 1 
both the longitudinal and transverse relaxation becomes 



3 



important, and it is difficult to find an analytical solu- 
tion to the FPE. Fortunately, in this temperature range 
cci and a 2 given by Eq. @ already become very close 
to each other, so that one can set ot\ => a 2 everywhere. 
Then the calculation yields 



i 



■ exp 



2 V a 



o-l 



xerfc 




(12) 



Here 9(x) is the step function and u>\ = jM e /(2xj_) is 
the ferromagnetic-resonance frequency. For a < 1 the 
total rate simplifies to 



r = 2a2(jJi 



exp 



-a 1- 



(13) 



This reduces to the Brown's formula T = 
2a 2 ^i •\/< J /7re _<T in the limit a — > 0. Exactly at 
a = 1, Eqs. {TTJ and (H2) yield V = a 2 LU 1 ^/aJire-' T/2 . 
Just below T c according to Eqs. (jBJ and © one has 
aa oc (T c -T) 2/3 ~ 7 , where /3 and 7 are the critical 
indices for the magnetization and susceptibility. Within 
the MFA aa remains finite at T Cl whereas for more 
realistic models it diverges. It makes, however, little 
sense to work out the appropriate limiting expressions 
for r because near T c the high-barrier approximation 
AU 3> T becomes invalid. In fact, the prefactor Tq 
in Eq. (jll|l does not strongly depend on a. The main 
difference of our result from the Brown's formula with 
a temperature-dependent barrier, AU oc M 2 (T) is 
described by the two factors: (i) additional lowering of 
the barrier because of the non-rigid magnetization, Eqs. 
(|10|l and H 1 1 p : (ii) crytical divergence of the damping at 
T c , Eq. 0. 

Brown has obtained the I /a correction to the escape 
rate for a ^s> 1 in the form T = 2a 2 LO\ \fo /Tre~ a (1 — 1/er) 
[lif . Within the LLB approach finding this correction in 
the whole range < a < oo is a complicated task. For 
a < 1 the correction factor in Eq. (|12Jl simplifies to 



(14) 



To illustrate the practical implication of our approach, 
we consider thermal switching of a FePt particle (mag- 
netic grain) at high temperature. The LLB-Langevin 
equation, Eq. has been integrated numerically with 
C = 0, (Ci) = 0, and 



2k B T 



«oM(t-0. (15) 



where i, j = 1,2 and /x, v = x, y, z. 

The microscopic relaxation parameter A in Eq. has 
been found analytically for a spin-phonon inter action 0- 
However it is difficult to obtain reliable theoretical results 
for A in general, as well as to extract A from experiments. 



For our illustration below we just set A = 0.1, neglecting 
its temperature dependence. The values of m e (T) of Eq. 
(@J can be measured or obtained from the Curie- Weiss 
equation m e = Bs(m e f3), where (3 = 5 2 Jo/(fcsT), Bs{x) 
is the Brillouin function, and Jo is related to the exper- 
imentally measured T c via T c = (1/3)5(5 + l)Jo within 
the MFA. For FePt T c = 750 K, and the best fit for 
m e (T) is obtained with 5 = 3/2 12]. For FePt we take 
M e (T = 0) = 1100 emu/cm 3 , K(T = 0) = 1.24 x 10 7 
erg/cm 3 , so that x± = = 0)/ [2K(T = 0)] = 

0.0488 Oe cm 3 /emu. In the same way, the longitudinal 
susceptibility \\\ = OM/dH can be measured or found 
analytically from the MFA: 



v M 2 e {T = 0) P&gimj}) 



5 2 J ( 



o I- PB' s (m e (3)' 



(16) 



where vq = 6.4 x 10 23 cm 3 is the unit-cell volume and 
B' s (x) = dB s {x)/dx. 

To integrate the LLB-Langevin equation, the Heun nu- 
merical scheme has been used. The physically rea- 
sonable interpretation of the stochastical process is that 
in the sense of Stratonovich, as was first stressed for the 
LL equation in Ref. 0|- Lately, it has been shown 
[l7j that even a naive Euler scheme which yields the Ito 
solution, would converge to the proper averaged phys- 
ical value, if the magnetization is normalized at each 
time step, reflecting the conservation of the magnetiza- 
tion length. However, in the case of LLB equation the 
magnetization length is also a stochastic fluctuating vari- 
able, so that the choice of the integration scheme should 
explicitly include the Stratonovich interpretation. 

The spins were prepared in the state n z = —1, and 
mean first-passage time (MFPT) time was evaluated as 
the time elapsed until the magnetization reached the 
boundary value beyond the barrier, n z — 0.5. The ex- 
act position of this boundary only slightly changes the 
MFPT. Alternatively, one can set the boundary at the 
top of the barrier, n z = 0. In this case one has to multi- 
ply the time by 2, since in 50% of all realizations the spin 
crosses the barrier and in 50% of all realizations it falls 
back [13]. In the high-barrier case, T < AU, the MFPT 
should coincide with the relaxation time T" 1 following 
from the FPE. 

For a comparison, we also solved the (naive) LL- 
Langevin equation with a constant but thermally reduced 
magnetization length, 

n = 7 [n x (H off + C)] - ia 2 [n x [n x (H cff + < 2 )]], (17) 

where H e g is given by Eq. J5J without the last term. 
The temperature dependence enters this equation, as in 
the LLB case, via the scaling of the anisotropy energy 
with M 2 (T) [Eq. and Eq. without the last term]. 
The non-rigorous derivation of Eq. I|17|) starts with the 
equation s = 7[s x (H cff + <)] - 7^[ s x [s x (H cff + £ 2 )]] 
for the spin vector of unit length s. [The same starting 
equation is used for the derivation of the LLB equation, 
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FIG. 1: Switching times for a FePt particle with 8nm diame- 
ter calculated numerically from the LL-Langevin and LLB- 
Langevin equations and analytically from the approproate 
Fokker-Planck equations. While the integral relaxation time 
(IRT) works very well at all temperatures, the high-barrier 
asymptotes break down for a < 1. The 1/<t corrections im- 
prove the asymptotes for a > 1. 



that is valid in the whole temperature range. The IRT 
for the LL model is also plotted in Fig. fI]showing a good 
agreement with the numerical data at all temperatures. 
The possibility to find the IRT analytically results from 
the spatial one-dimentionality of the FPE in the axially- 
symmetric case: / = f(9, t). For the LLB model there are 

two spatial coordinates, {9, n} or |n 2 , n± — ^Jn^. + riy j , 

and a rigorous analytical solution for the IRT seems to 
be impossible. 

In conclusion we have introduced the formalism of 
the temperature fluctuations within the mean field ap- 
proach based on the Landau-Lifshitz-Bloch equation. 
The new Langevin equation could serve as a basis for 
the temperature-dependent micromagnetic approach for 
small particles (or discretization elements) at high tem- 
perature, similar to standard temperature-dependent mi- 
cromagnetics but valid in the whole temperature range. 
This new micromagnetics will have practical importance 
for the heat-assisted magnetic recording applications. 



Eq. JQ), in the classical case]. Replacing s in this equa- 
tion by its thermal average, s m = (s) , and rescaling 
m as m = nm e yields Eq. (|17fl with a 2 — \m e . The lat- 
ter is in a striking contradiction with the rigorous LLB 
expressions for the damping parameters, Eq. |J3J. This 
difference becomes very important at elevated tempera- 
tures, and there is no easy way to improve the naive LL 
results. In our simulations of Eq. (|17|) we just use the 
constant A = 0.1 instead of qg, to conform with existing 
publications, e.g., with Ref. .12]. Using a2 = Am e leads 
to even more pronounced difference between the LL and 
LLB results. 

Fig. fl] shows the MFPT of a 8 nm one-domain FePt 
particle as a function of temperature, calculated numer- 
ically from the LL-Langevin and LLB-Langevin equa- 
tions. These numerical results are compared in Fig. Q] 
with Brown's analytical formula for the relaxation time 
T _1 0] and the result of Eq. lJT2)l. The switching time 
calculated within the LLB approach is always lower than 
that of the LL approach due to the additional lowering 
of the energy barrier (|10fl and the critical growth of the 
damping at T c . We have also shown the temperature 
T* ~ 738 K at which a = 1 and the geometry of the 
barrier changes. For a given particle's volume, our high- 
barrier approximation leading to Eq. I|llfl becomes in- 
valid for T > T*. We cannot increase the volume, how- 
ever, without violating the single-domain requirement. 

Both Brown's formula for the LL model and Eq. 111|) 
for the LLB model describing the Arrhcnius dynamics 
are only valid for T <C AU. Switching times showing an 
unphysical divergence near T c is the signature of their 
breakdown. For the LL model, there is a better analyt- 
ical approach describing the thermally activated escape 
in terms of the integral relaxation time (IRT) 0, 0, E| 
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